Spectral Effects of Strong Chi-2 Non-Linearity for Quantum Processing 
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Optical x' 2 ' non-linearity can be used for parametric amplification and producing down-converted 
entangled photon pairs that have broad applications. It is known that weak non-linear media exhibit 
dispersion and produce a frequency response. It is therefore of interest to know how spectral effects 
of a strong x' 2 ' crystal affect the performance. Here we model the spectral effects of the dispersion 
of a strong x' 2 ' crystal and illustrate how this affects its ability to perform Bell measurements and 
influence the performance of a quantum gates that employ such a Bell measurement. We show that a 
Dyson series expansion of the unitary operator is necessary in general, leading to unwanted spectral 
entanglement. We identify a limiting situation employing periodic poling, in which a Taylor series 
expansion is a good approximation and this entanglement can be removed. 

PACS numbers: 03.67.Lx, 42.50.-p 
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I. INTRODUCTION 

An optical x^ non-linearity can combine two lower 
energy photons into one higher energy photon via the 
parametric up-conversion process, and conversely, break 
a higher energy photon into two lower energy photons 
via parametric down-conversion, such that the total en- 
ergy of the photons is conserved before and after the 
conversion. Optical x non-linearity is widely used for 
parametric down conversion experiments to produce po- 
larization entangled photon pairs, as well as for paramet- 
ric amplification. Entangled photons have a large variety 
of applications, such as demonstration of Bell's inequal- 
ities violation quantum error encoding produc- 
tion of heralded single photon sources 9] , quantum tele- 
portation quantum dense coding 5[, and entangle- 
ment swapping Q . In addition, in principle a sufficiently 
strong x non-linearity could be used to perform deter- 
ministic Bell measurement 0], as we shall discuss later 
in this paper. Bell measurements can be used in appli- 
cations such as, teleporting qubits 0], transferring quan- 
tum information with quantum repeaters as well as 
performing quantum computation [9(. 

The conversion efficiency of current x^ non-linearity 
is far below unity and past research has mostly concen- 
trated on the properties of weak x ■ However, the 
strength of x' 2 ' has been improving, as demonstrated 
in high photon number experiments, and it is therefore 
of increasing interest to examine the properties of strong 
X*- 2 ) media. Furthermore, since x' 2 ' media have intrinsic 
spectral response, it is important to understand how the 
spectral effects of a x^ medium affects the rate of up and 
down conversion. S pec tral effects of weak x' 2 ' crystals 
have been examined [lfj and some numerical research has 



been done for x^ non-linearity of arbitrary strength 11 1 . 
In this paper, we shall model the spectral properties of a 



strong x'- 2 -' crystal and explore limits where analytical so- 
lutions can be found. To do so, we shall examine how the 
dispersion of a strong \ non-linear crystal affects the 
profiles of spectrally Gaussian input photons and deter- 
mine how the probability of up-conversion depends on the 
dispersion. Specifically, here we consider Type II conver- 
sion. Moreover, we put forth a Bell measurement scheme 
based on parametric up-conversion and further develop 
a quantum gate from it, and examine the success rate 
of the gate under the influence of the dispersion of the 
crystal. Note that although our discussion focuses on up- 
conversion, a strong x' 2 ' non-linearity has both up and 
down-conversion happening at the same time and thus 
the conditions and results that we derive here would also 
similarly apply to down-conversion. 

This paper is arranged in the following way. The next 
section discusses how we model the spectral properties of 
a x' 2 ' non-linear crystal. We indicate the problem of non- 
commutivity of the interaction Hamiltonian at different 
times, which requires us to use the Dyson series 12, Hf to 
calculate higher order effects of the crystal on the evolu- 
tion of the photon states, as opposed to using the simpler 
Taylor series. The Dyson series leads to spectrally mixed 
states and this is undesirable for most applications. On 
the other hand, as we will show, spectrally separable solu- 
tions exist for the Taylor series. In section III, we exam- 
ine the case where we send in a pair of separable photons 
through slices of weak \^ crystal that are well sepa- 
rated and derive the output state for the photons and 
the probability of up-conversion. Then in section IV, we 
examine the case where we have periodically poled bire- 
fringent x^ spacers in the crystal, and likewise derive 
the output state and the up-conversion rate. A rather 
surprising result of this paper is that we predict the ef- 
ficiency of a strong bulk x^ crystal, where the unitary 
evolution is modeled using the Dyson series, will be lower 
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than the efficiency of the aforementioned thin slices case 
and the periodical poling case of a strong x*- 2 -* medium, 
where the unitary evolution is modeled using the Taylor 
series. In section V, we describe how one may construct 
a quantum gate based on Bell measurement with strong 
non-linearity and find the probability of success of 
the gate. We conclude in section VI. 



II. MODELLING THE SPECTRAL PROPERTIES 
OF A CHI-2 NON-LINEAR MEDIUM 

The Hamiltonian of a non-linear x^ medium is spec- 
trally dependent and it is crucial to understand how the 
spectral response of the medium affects the conversion 
rate, as well as the spectral profile of the output photons. 
So in this section, we shall model the spectral properties 
of a x^ medium. The process of x non-linearity can 
be studied in the interaction picture and the unitary evo- 
lution of a state vector is given by: 



W) (i) 



where T is the time ordering operator. Taking into ac- 
count the time ordering leads to the following Dyson se- 



ries 
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If the interaction Hamiltonian commutes at different 
times, then the time ordering operator in equation (TTJ) 
has no effect and can be dropped, resulting in the usual 
Taylor series for the unitary expansion. The interaction 
Hamiltonian for a x^ process has the form: 



nj(u!j) is the refractive 



where Aj(u>j) = i. 

index for mode j and S is the cross section area of the 
beam. We assume that Aj(uij) = Aj is slowly varying 
for the frequencies of interest, allowing it to be factored 
outside the integral. The frequency integrals have lower 
bounds extended from zero to negative infinity. This is 
mathematically legitimate because we are considering a 
system that operates at high frequency, where essentially 
there is no population present at low frequency. When 
we further integrate the Hamiltonian over z from to L, 
we obtain: 



H(t) =xL JJJ di^ p du s duia) p (uj p )d s {u} s )di(uji) 
. ( LAk\ 
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where Ak = k p {oj p ) — k s (u) s ) — fcj(wj) is the phase mis- 
match and Aw = lu p — lu s — u>i is the frequency detun- 
ing. Here x is again the interaction strength but incor- 
porated with some constants from the electric field ex- 
pressions. Following Grice and Walmslcy [HI, we Taylor 
expand the phase mismatch and retain terms up to first 
order by assuming higher order terms arc negligible, then 
Afc ?s AA;*- ' + k'pi/p — k' g v s — k'^i, where i/j = toj — fj,j and 
fij is the centre frequency of the photon in mode j. We 
set /i s = /Zj = n and \x v = 2/i. The parameter k'j is 
the derivative of wavenumber kj with respect to uij and 
evaluated at Hj . Due to conservation of momentum, the 
zeroth order term Afc(°) = k p (fx p ) — fc s (/t s ) — ki{jJbi) = 
and thus Ak « k' p v p — k' s v s — k[vi. In our calculation, 
we assume that k' s — k' i ^ 0, which is the case of Type II 
parametric conversion. 

It can be shown that the interaction Hamiltonian in 
equation (|5|) docs not commute at different times and 
therefore the Dyson series applies, instead of the Taylor 
series, when calculating higher order terms in the unitary 
expansion. The first order terms are identical, so we shall 
quantify and compare the similarity between the second 
term of the Taylor series and the Dyson series. A single 
photon with a Gaussian spectral profile is given by: 



t 



H{t)=x w e dv 6 E p {r,t)E a {v,t)Ei{v,t)+h.c. (3) 
Jv 

Hereafter, we shall simplify the analysis to one spatial 
dimension, the propagation direction. This is legitimate 
if we consider collinear Type II conversion or in the case 
that the net transversal effects arc negligible. The ex- 
pression for the electric field operator of mode j with 
spatial degree z is: 



du)jAj(LOj)df(ujj) exp(ikj(u>j)z — ujjt) 

(4) 



dwj d\u>j)f(ujj)\0} 



(0) 



where f{u)j) = exp exp (if i/,-) and 

(1|1) = / \ f(ujj)\ 2 dujj = 1. Hence for up-conversion, the 
spectrally separable two-singlc-photon input state with 
modes s and i is: 



dwgdcji a\ (oj s )al(u>i)f(ui s )f(uji)\0) (7) 



Let \ip2,T) and \ip2,D) be respectively the state compo- 
nents arising from the second term of the Taylor series 



3 



and Dyson series expansions of the unitary operator act- 
ing upon the input state. Mathematically, 
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and 



_ (fc'j-fcl) (fc> s +fc>i) / k' z i +k' i 3 k' s +k' i k' s 3 +k[ 4 

9 ~ 2V2(k'3+k' s 3) y CT 2 (fc^-fe^) 2 ( fe ; 2 + ^ 2 )- 
To quantify the similarity of the states IV^.t) and \ip2,D), 
we define their fidelity to be: 



F, = 



2.D 



VI(V'2,t|V'2,t)||(V'2 : d|V'2,_D>| 



(10) 



For the set of reasonable parameters, k' s = 5.6 x 
10- 9 (s/m), fc£ = 5.2 x 10~ 9 (s/m), and cr = 
10 9 (i/z), and assuming the extended phase match- 
ing condition k' p = i , and the special condition 
L 2 7 cr 2 (k' s - k'p) (k' p - kfi = i, we obtain F 2 = 0.747, 
which means the second term of the Dyson series sig- 
nificantly differs from the second term of the Taylor se- 
ries. From equations ([5]) and ©, it is clear that this 
difference comes from the Erf function in \i})2.d)- By as- 
suming the extended phase matching condition and the 
special condition, |"02,t) can become spectrally separa- 
ble [13] and proportional to \ipo), however, the Erf func- 
tion induces spectral entanglement between the two pho- 
tons of IV^.d), which makes the two states substantially 
different. In principle, it is possible to test this difference 
experimentally by examining the efficiency difference be- 
tween a moderately strong bulk x crystal and for the 
same crystal being cut into many thin slices and sep- 
arated sufficiently apart. We shall discuss the latter in 
the next section. 



From the complexity of the 2nd order term of the 
Dyson series, we doubt that the high order terms in the 
series can have the spectral entanglement canceled out. 
Besides, for terms higher than the first order, the upper- 
bound of the time integrals is a time variable that has to 
be integrated by the next time integral. This makes the 
calculation of the sum of the Dyson series very difficult 
and complicated. If x^ is sufficiently weak, then the 
higher order terms may be neglected and the calculation 
is tractable in this limit. For instance, Grice and Walms- 
ley (Io| have examined weak parametric down conversion 
by ignoring higher order terms in the series. However, 
here we are interested the spectral effects of a strong x^ 
crystal and we must include higher order terms for cal- 
culating the evolution of the input states. In the next 
section we shall look at the case where the strong x^ 
crystal is cut into thin slices and separated sufficiently 
apart. Since the interaction is weak, we may legitimately 
ignore higher order terms for the evolution of the state 
provided by each thin slice. 



III. OBTAINING STRONG CHI-2 
NON-LINEARITY FROM MANY WEAK SLICES 

In this section, we examine the case where a strong x 
non-linearity is obtained from many thin slices separated 
sufficiently apart, such that the wavepacket exits one slice 
before entering another and each slice provides only a 
weak interaction. In practice, this may be accomplished 
by having the photons in each mode passing one thin slice 
for many times in a loop. If the total length of a bulk 
X^ crystal is NL and we divide it into TV pieces of equal 
length L, then the unitary operator can be re-expressed 
as: 
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H(t)dt 



-AT 



H(t)dt 



H(t)dt 



(11) 



where AT = (t\ — to)/N. Each factor in the expression 
represents a weak interaction by a slice of crystal. There- 
fore, if the slices are sufficiently apart, the time bounds of 
the integrals can be conveniently extended to — oo and oo 
by considering far field limits [lj], and thus each factor 
in the expression is the same. The unitary operator can 
now be expressed as a Taylor series and the time ordering 
operator becomes irrelevant: 
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N 



cxp ( — / H{t)dt 



(12) 



By integrating the interaction Hamiltonian in equa- 
tion over time t, we have: 

/>oo 

H = H{t)dt 



=XL JJJ duj p du s duj l dl(uj p )a s (uj s )di(uji) 
f LAk\ fiLAk' 



sine 



T 



cxp 



S(Atu) + h.c. 



=H+ + H- 
where we define 



H+= X L 



(13) 



dwpdujs dujt a P (ujp)a s (uj s )ai(uJi) 



(14) 



H- =X*L JJJ dui p duj s dui l ap(u; p )dl(uj s )al(uji) 

§*(u a ,Ui)5(Au) (15) 



fLAk\ fiLAk\ , . 
$(w s , uji) =smc I — ^— j cxp I — — j (16) 



Using equations (JT2J) to (jT6J), we shall now derive the 
evolution of the input state \ipo}- To calculate the first 
order term of the Taylor series, we act the Hamiltonian 
from equation (|13| onto the input state, which gives: 



^_H\^ ) = ^H + \ip ) = ^f <kj p al(u p )J p \0) (17) 

where J p = / duj s f(uj s )f(ujp ~ lo s )$(lu s , lu p - uj s ). Sim- 
ilarly, the second order term of the Taylor series is: 



4® 



1 fN 



2! \ih 



_ f X L ^ [ [ dLd s duj i a i s (uj s )dl(Ld l )§*(uj s ,uJi)J s s\0) 



2! V ih J JJ 



(18) 



where we further define: 

Likewise the third order term of the Taylor series is: 



1 /AT 
3! \ih r 

3! I ih ) 



du) p a P {ujp)R p Jp\ti} 



(19) 



where R v 



J duj\§(u 



u>)\ 2 . Each of these terms 



contains the function </, which is obtained by integrat- 
ing the product of the Gaussian spectral profiles of the 
photons and the sine frequency response function of the 
%( 2 ) medium. Since it is difficult to integrate prod- 
ucts of Gaussians and sine functions, in our calcula- 
tion, we have made use of the approximation, sinc(x) « 
y/jTr cxp(— 7X 2 ), where the parameter 7 m 0.193..., is 
derived from equating the full-width-half-maximum of 
the two functions. After integrating the expression of 

R p , we get R p = / — f~~T^i which is independent of 

y L2 { k 's~K) 

Lj p . Thus we can drop the p subscript and set R = R p . 

Lets suppose \ip e ven) = cos(NH/h)\ip ) and \ip dd) = 
s'm(N H / h)\ipo) are, respectively, the sum of the even and 
odd terms of the unitary expansion after acting on the in- 
put state, such that the output state is \ipout) — U\ipo) — 
iV'euen) - * I VWd) • The odd state represents the part of the 
output state in which the two photons in the signal and 
idler modes are up-converted into a photon in the pump 
mode. The even state represents the part of output state 
in which the two photons in the signal and idler modes 
are not up-converted and remain in the two modes. From 
the definition of the odd state, we obtain: 



\ipodd) =sin ( —H ) \ip ) 



jo 



B 



du> p a P (ujp) —¥= sin 
v R 

X\NL 



X \NL 



R)\0) 



R 



sin 



R \^p) 



(20) 



where B = 



2+L 2 ~ l a 2 (k' s -k' i ) 



J1 P fc A2 , H'p) = J duj p dl(ujp)fp\0}, 



■fp = J pk„ ex P(-I^)> 9 is the argument of 
the complex number x m polar form, and a p = 

j a 2 (2+L 2 7 o- 2 (/^-/cQ 2 ) 

V l + L2 ~t° 2 ({K-K) 2 +(K-KY) 
of the even state gives: 



Similarly, the definition 
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\lpeven) = COS ( ) \lp ) 



1 (\ X \NL 



dtjj s duJia\(L0 s )a\ {u>i ) 



$*(w s ,cj i )J s , i + ---J|0) 



Note that if fc- -> k' then i? 



277T3 



(21) 



CO, 



i2 ( fc L-M) 

and the sin function in the odd state will have an un- 
physical infinite oscillation. This is caused by ignoring 
higher order terms of Ak. When k[ — > k' s , the contribu- 
tion from higher order terms of Ak becomes more signif- 
icant and the infinity will be prevented if these terms arc 
included. Nevertheless, in the following, wc shall show 
that the probability goes to zero when k[ — > k' s , which 
means Type II conversion simply does not work in that 
limit and thus there is no observable physical problem. 

With the expression for \ipodd), we can now determine 
the probability of having a pump photon by calculating 



P(odd) 



l/ (lpodd\lpodd) 



\(4>odd\4>odd}\, which gives: 



B 2 
R 



(22) 



Setting l -Z^±VR = § gives: 



P{odd) 



R 



4(d s -diY 



^ (l + d2 + d?)(2 + (d a -di) 2 ) 



(23) 



where d s — L^fa {k' s — k' p ) and di — L^/ja (fc£ — k' p ) . 
Figure [T] shows the plot of P(odd) against d s and di. The 
maximum probability is one and occurs at two points, 
di = —d s = which is achieved when wc have the 



, as well as 
The 



extended phase matching condition k' p = ^ * 
the special condition L 2 ja 2 (k' s — k' p ) (k' p — k^) = |. 
plot also shows that the probability is zero at di = d s , 
that is when k' s = k[. Beware that the zero probability 
trough is only meaningful in the sense that k' s is close to 
k[ but not equal. This is because, strictly speaking, in 
order to understand what happens at k' s = k'^ we have 
to include higher order terms of Ak when calculating the 
probability. 

At the optimal points, the extended phase matching 
condition and the special condition are satisfied, and 
as we mentioned in the previous section, these condi- 
tions allow the second term of the Taylor series to be 
spectrally separable and proportional to IV'o)- In fact, 
these two conditions make ^*(ui 8 ,U)i)J Si i = Rf(u) s )f(u>i), 
and thus all even order terms in the series are spec- 
trally separable and proportional to \ipo), such that 



\i>even) = cos ( \/RJ \ipo) ■ From equation (|20p . we 
can see that all the odd terms are proportional to \tp p ) as 
iV-Wd) = e' 9 sin ^ ^ X ^ L VRj IV'p)- Physically this means 
that the evolution is a Rabi oscillation between the only 
two possible basis states \ip p ) and IV'o)- So inside the 
crystal, two processes happen concurrently, the two pho- 
tons are up-converted into a pump photon and the pump 
photon is down-converted exactly back to the original 
two photons. Hence, if we choose the extended phase 
matching condition and the special condition, spectral 
entanglement between photons is avoided, and moreover, 
no population can be leaked to other spectral modes. So 
by tuning the Rabi oscillation such that ^ L VR = f, 
deterministic up-convcrsion can be achieved. 



P (Odd 



0.75 




FIG. 1: (Color online) Plot of the success rate in converting 
two photons, one in the signal mode and one in the idler 
mode, into a single photon in the pump mode using strong 
Type II parametric up-conversion, where V^R = \- The 

rate is plotted against dimensionless parameters di and d s . 



The probability of successful conversion will be low- 
ered if the extended phase matching condition and 
the special condition are not perfectly achieved. Lets 
suppose that some errors, t\ and £2, are allowed in 



the conditions, such that k' p = 

(i + £2 ) 2 iV 2 (K-K) iK-K) 

\X\NL 



(1 



1 2 

Assuming that 



R 



-c — y — ^ , and having the set of reasonable param- 
eters, k' s = 5.6 x l(T 9 (s/m), k'i = 5.2 x lCT 9 (s/m), and 
(7 = 10 9 {Hz), then for e x = 0.01 and e 2 = 0.01, the 
probability of success is 0.9803; and for ei = 0.001 and 
62 = 0.001, the probability of success is 0.9998. Hence 
high conversion rates can be achieved if the errors in the 
extended phase matching condition and the special con- 
dition are reasonably small. 
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IV. STRONG CHI-2 NON-LINEARITY MEDIUM 
WITH PERIODIC POLING 

In the previous section, we analyzed the spectral ef- 
fects of a x medium obtained from many thin slices of 
crystal, such that the contributed non-linearity of each 
slice is small. This allows us to apply the simpler Taylor 
series to the unitary operator expansion, instead of the 
complicated Dyson series. However, guiding the modes 
through many thin slices, or many times through a sin- 
gle slice, could be experimentally challenging. In this 
section, we examine an alternative situation where we 
have a bulk medium comprising TV slices of x^ crystal 
and with TV — 1 spacers in between. We prove that in the 
limit of sufficiently large TV, it is valid to use the Taylor 
series to approximate the Dyson series. Figure [5] shows 
the schematic of TV slabs of x' 2 ' crystal of length L and 
TV — 1 slabs of x' 1 ' spacers of length h in between. 



• • • 




h h 



h 



FIG. 2: Schematic of TV slices of x^ crystal of length L, 
with TV — 1 x spacers of length h in between. 



U'Ren et al [H have shown that the interaction Hamil- 
tonian of a medium consisting of TV x^ crystals and 
TV — 1 periodically poled x spacers has the expression: 



H(t) = X L 



sine 



/ LAk\ iL&k sin(^) nN-m 



e 2 



sin(f) 



-e 2 



h.c. (24) 



where <f> = LAk + hAn and An = K v v v — K s v a — KiVi is the 
phase mismatch introduced by each of the birefringent 
spacers, where Hj is the wavenumber for field j, taking 
into account the dispersion of the spacers. Rewriting the 
sin functions as sine functions and approximating them 
as Gaussian functions gives: 



H(t) = X NL 



dLO p du s dio i al(ujp)a s (uj s )a. l (uj. l )t 



fLAk\ 



-LAfc -1 a2 ; (jV-jj 



V 2 J 



l + + h. 



— i/S.Lot 



(25) 



If hn p = Lk' p , tiK s = Lk[ and /ire, = Lk' s , then 
4> = L(2k ' v p - (k' s + k'^ivs + Vi)). Further applying 



the extended phase matching condition gives </> = r/Auj, 
where rj = L(k' s + k[). Now the interaction Hamiltonian 
becomes: 



H(t) = X NL 



du)pdw s du>i aUujp)a s (uj s )ai(u)i)<: 



sine I — — | e 2 e 7 4 ' e 2 ^ + h.c. 



(26) 



Recall that the interaction Hamiltonian in equation (fS]) 
does not commute at different times, which has to do 
with the sine function in the expression. In the limits 
r/\/ N 2 — 1 ^> Lk'p where j is for all three of the modes, 
the sine function in equation (|26p is relatively flat in the 
domain where the Gaussian function is significant. This 
means that sine (^j^) ~ 1 within the domain. Further- 
more, the phase exp ( lL ^ k ) can be neglected in this limit. 
Hence the interaction Hamiltonian tends to: 



H(t) « X TVL 



dw p du! s duj l a)p{ijjp)a s (uj s )a l (uji)(. 



Surprisingly, this interaction Hamiltonian commutes 
at different times and we may again use the Taylor se- 
ries to calculate higher order terms of the unitary evo- 
lution. Using the extended phase matching condition, 
one of the three limits, \/ N 2 — lrj ^> Lk' p , implies that 

TV > ^ w 1, thus we expect that for TV > 1 the Taylor 
series should give a good approximation to the Dyson se- 
ries. To confirm this and determine how large TV has to 
be in practice such that the Taylor series gives a good ap- 
proximation to the Dyson series, for various values of TV, 
we calculate the parameter F2 that we defined in equa- 
tion pTT)) . Figure [3] is the plot of F2 against TV. It shows 
that the fidelity between the 2nd order term of the Taylor 
series and the Dyson series is 0.998 for TV = 5 and con- 
tinues to increase as TV gets larger. We have also checked 
that the phase difference between the two terms is negli- 
gible. Hence we argue that for a sufficiently large number 
of spacers placed between the x^ crystals, we may use 
the Taylor scries to approximate the Dyson series. In- 
tegrating the interaction Hamiltonian in equation (|26p 
over time and simplifying the equation, we can write the 
Hamiltonian as: 



H =xNL 



dwpdu} s dwia p ((jjp)as{u) s )ai(uji) 



( LAk\ iLAk 
smc I I e 2 d(Auj) + h.c. 

--NH 



(28) 
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Since the Hamiltonian is for all slices of crystals and 
spacers, the state \ipodd) = sin(7Y/^.)|V , o)- Hence, 
the profile of the up-converted photon is \tp dd) = 
e %B -!j= sin ( ^ X ^ L V~Rj and the probability of re- 
conversion is P(odd) = -^r sin 2 O x ^ L \/~Rj, which are 

the same as the profile and probability that we have 
found in the case of many thin slices. So if we choose 
the extended phase matching condition and the special 
condition, which maintain spectrally separable photons 
and gives Rabi oscillation between \ip p ) and |V>o) basis 
states, then by tuning the oscillation to ^ X ^ L >/R = -|, 
deterministic up-conversion can again be achieved. 



| Target) - 
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FIG. 4: Schematic of the Gottesman and Chuang 
teleportation CNOT gate. The resource state 

\Resource) = {{\HH) + \VV))\HH) + (\HV) + \VH))\VV))/2 
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FIG. 3: Plot of F 2 against N 



V. BELL MEASUREMENT AND QUANTUM 
GATE SCHEMES WITH STRONG CHI-2 
NON-LINEARITY 

Gottesman and Chuang 0] showed that it is possible 
to build a CNOT gate by means of quantum telepor- 
tation and post-selection. Figure [4] shows the CNOT 
gate scheme. The control qubit \control) = a\H) +b\V) 
and the target qubit \target) = c\H) + d\V) are 
two arbitrary single qubit states in polarization en- 
coding. The resource state is the entangled state 
((|M> + \VV))\HH) + (\HV) + \VH))\VV))/2, which 
is prepared offline. The two input qubits are each sent 
to a separate device to perform a Bell measurement with 
the resource state. After measuring the input qubits 
and two of the four qubits of the resource state, the 
results are then used to perform single qubit operations 
on the two remaining output qubits. These single 
qubit operations can be done with optical waveplates. 
Gottesman and Chuang showed that this procedure 
teleports the input state to the output qubits with a 
CNOT operation applied between them. 

The problem with implementing this gate with linear 
optics is that the Bell measurements only work 50% of the 



time, giving a | probability gate. Boosting the probabil- 
ity is possible but involves complex circuits [l5| . To per- 
form the Bell measurement deterministically, we propose 
to use a sufficiently strong x^ crystal along with linear 
optics as shown in FigureO The strong x material con- 
sists of two slabs. One slab performs a Type II paramet- 
ric up conversion for \H) s \V)i and the other for \ V) s \H)i. 
One may also construct a Bell measurement circuit us- 
ing non-collinear Type I parametric up conversion in a 
similar fashion. The purpose of a Bell measurement is to 
distinguish the four different Bell states. To see how the 
scheme works, lets consider the cases of inputting the four 
Bell states. The two Bell states (\H) s \V). l ±\V) s \H). l )/V2 
are up-converted at the x^ crystal to a pump photon 
with (\H) P ± \V) p )/y2 states respectively by each slab. 
Thus measuring the polarization of the pump photon al- 
lows us to distinguish the two Bell states. If the conver- 
sion fails, the two photons will pass through the crystal 
and the half wave plate (HWP) in the i mode will turn 
the states into (\H) s \H) l ± \V) s \V)i) / ^2. These two 
states will have the photons bunched at either one of the 
two 50:50 beam splitters, giving two photons arriving at 
either one of the four detectors, D\, D2, -D3 or D4, and 
thus signaling a failure event. Regarding the other two 
Bell states (\H) s \H)i ± \V) S \V)i)V2, they simply pass 
through the x crystal and the half wave plate in the i 
mode will turn the state into {\H) s \V) i ±\V) s \H) l )^/2. 
The photons will then anti-bunch at the two 50:50 beam 
splitters. For {\H) s \V) l + \V) s \H)i)y/2, the two photons 
will arrive at either D\ and D3 or at D2 and D4. For 
(\H) s \V)i — \V) s \H)i)y2, the two photons will arrive at 
cither D\ and D4 or at D2 and D3. Thus, identifying the 
combination of flagged detectors allows us to distinguish 
the remaining two Bell states [i~8j . 

Although each Bell measurement uses two slabs of x^ 
crystal, the probability of conversion is P(odd) and not 
P{odd) 2 because there are only two photons undergo- 
ing one up-conversion process at any one time. Since 
two of the four Bell states are identified with a prob- 
ability of P(odd) and the other two arc identified de- 
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tcrministically, the probability of a successful Bell mea- 
surement is therefore | ^1 + sin 2 y x ^ L VRj J • The 
CNOT gate comprises two Bell measurements that both 
need to succeed, so the probability of success of the gate 

is P 2 (odd) = \ (l + ^sin 2 (kL^Vi?)) 2 . Given the 
extended phase matching condition and the special con- 
dition, and setting ^ X ^ L VR = § , the success rate of the 
Bell measurement and the gate becomes one. Hence by 
employing sufficiently strong x' 2 ' crystals with the cor- 
rect phase matching conditions, Bell measurements can 
be done with unit probability, and a deterministic CNOT 
gate can be achieved via teleportation. We note that Kim 
et al has proposed a similar Bell measurement scheme 
that requires two x crystals and four slabs in total, to 
perform a Bell measurement. Since our scheme requires 
one crystal less, it gives a higher success probability over 
Kim et al's scheme if the conversion efficiency is less than 
one. 




FIG. 5: Schematic of the deterministic Bell measurement 
setup. The \^ medium does TYPE II up conversion. 
The three modes s,p, i are shown as different non-collinear 
spatial modes for the purpose of illustration. However, this 
is equivalent to collinear TYPE II conversion if one makes 
use of polarization beamsplitters and dichromatic beamsplit- 
ters to combine the three spatial modes into one before the 
conversion and then separate them again after the conversion. 



IV. CONCLUSION 

In this paper, we have modeled the spectral effects of 
the dispersion of a x' 2 ' crystal and pointed out that be- 
cause the interaction Hamiltonian does not commute at 
different times, higher order terms of the unitary evolu- 
tion of the state ought to be calculated using the Dyson 
series instead of the Taylor series. We have quantified the 
similarity between the second order term component of 
the state for the two series and estimated that the fidelity 
is i*2 = 0.747. Also, we have shown that the Dyson se- 
ries induces spectral entanglement between the two input 
photons. We argued that this indicates that the conver- 
sion efficiency of a bulk x^ 2 ' crystal will be low. If the 
crystal is cut into many thin slices and well separated 
apart or has spacers between the slices, then the unitary 
operator can be calculated using the Taylor series. The 



two input photons can remain spectrally separable for the 
Taylor series, and thus high efficiency can be obtained. 
For the case of having N — 1 x spacers periodically 
poled between A slices of \^ crystals, we argued that 
the interaction Hamiltonian approximately commutes at 
different times in the limit of large A~ and demonstrated 
how the Taylor series is a good approximation to Dyson 
series in the limit by calculating Fi for various values of 
N. We have estimated that with N = 5, the fidelity F2 
is roughly 0.998 and the fidelity continues to improve as 
N increases. We have derived the up-conversion rate to 
be P{odd) = ^ sin 2 (^^Vr) for both the many thin 
slices case and the spacers case. We found that if the 
extended phase matching condition k' p = fc "^ fci and the 
special condition L 2 -fcr 2 (k' s — k' p )(k' p — k^) = | are sat- 
isfied, and that ^ L VR = § , then the conversion rate 
becomes unity. We have further discussed how a Bell 
measurement can be done by employing a strong x 
non-linearity to perform up-conversion on polarization 
encoded photons, and further implement it for construct- 
ing a teleportation-type CNOT gate. We calculated that 
the probability of success of the Bell measurement is 
i ^1 + sin 2 ^ ^fr L VR^j ^ , and the success rate of the 

CNOT gate is P 2 (odd) = \ (l + ^sin 2 (kL^^)) 2 . 
Again, by having the extended phase matching condi- 
tion and the special condition, and that ^f* L VR = f, 
for both the many thin slices case and the spacers case, 
deterministic Bell measurement and CNOT gate can be 
constructed. In this paper, our discussion has focused on 
up-conversion, however a strong x^ non-linearity has 
both up and down-conversion happening at the same 
time and thus the conditions and results that we derive 
here would also similarly apply to down-conversion. 
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